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Abstract 



We review the time evolution of wavepackets at the metal-insulator tran- 
sition in two- and three-dimensional disordered systems. The importance of 
scale invariance and multifractal eigenfunction fluctuations is stressed. The 
implications of the frequency- and wavevector-dependence of the diffusion co- 
efficient are compared with the results of numerical simulations. We argue 
that network models are particularly suited for the investigation of the dy- 
namics of disordered systems. 
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I. INTRODUCTION 



The time evolution of wavepackets reflects clearly the localization properties of the 
eigenstates of a system. While in infinite metallic systems wavepackets spread indefinitely 
and take on a Gaussian shape for long times, diffusion is absent in localized systems and 
wavepackets stay in a finite region around their starting point. In fact, the consideration 
of the diffusion of wavepackets has been the starting point of the theory of localization in 
random media JT|. At the mobility edge between metallic and insulating behavior the time 
evolution of wavepackets will be intermediate between the two extreme cases discussed above. 
The dynamics of wavepackets can be described by a frequency- and wavevector-dependent 
diffusion coefficient D(q,u). Its behavior has previously been studied in detail |g-f!|. The 
diffusion coefficient exhibits two characteristic features. First, from the scale invariance of 
the conductance || follows that the conductivity and via the Einstein relation the diffusion 
coefficient are scale dependent in higher than two dimensions. Thus only in two dimensions 
there exists a finite limit of D(q,u) as frequency and wavevector tend to zero. Secondly, on 
short distances and frequencies the diffusion coefficient reflects the multifractal correlations 
within single eigenfunctions of the system ||. 

In this paper we focus on the time evolution of wavepackets. While the same informa- 
tion is encoded in the frequency- and wavepacket-dependence of the diffusion coefficient, 
we concentrate on the time dependence of characteristics of wavepackets like their return 
probability, their moments, and their shape. Our purpose in doing so is twofold. On the 
one hand elucidate these discussions some of the crucial differences between time evolution 
in the diffusive regime and at the mobility edge. On the other hand can the time evolution 
of wavepackets directly been studied in numerical simulations of network models, allowing 
for a check of predictions of scaling theory. 

The paper is organized as follows: we first summarize the scaling form of the diffusion 
coefficient D(q, u). Next we present the time dependence of the moments, return probability, 
and shape of wavepackets that results from this scaling form. And finally, we compare these 
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predictions with numerical results obtain for network models. 



II. SCALING FORM OF THE DIFFUSION COEFFICIENT 

In this section we will briefly review known results about the scaling form of the diffusion 
coefficient. For a more complete discussion see, e.g., refs. 

Consider the time evolution of a wavepacket that at time t = is localized as a 5-function 
at the origin. Due to conservation of the probability the Fourier (Laplace) transform of its 
probability density with respect to space and time has of the form 

— iu + D(q, ijj)q l 

A scaling form for D(q,cj) that is consistent the conductance being scale invariant at the 
metal-insulator transition is 

D(q,u) = ^-F(Z/L u ,Zq), (2) 

where £ is the localization length, p = p(E) the density of states near the critical energy 
E = E c , F is a scaling function, and 

L u := { P nu)- 1 ' d = Lifnu/A)- 1 ^, (3) 

is a third length scale, besides g _1 and £, relevant at the critical point (A = (L d p) _1 is 
the mean level spacing). At the critical point, where the localization length £ diverges, the 
diffusion coefficient can be written as 

D(q,u>) = £lf(qL u ). (4) 

From various limiting cases the following form of f(x) can be deduced: 

c a ■ x {2 - d)/d : x -> 

cp : for intermediate values of x (5) 
c 7 • x~ v ^ d : x —> oo 



/(*) 



The behavior for large x is governed by multifractal density fluctuations in single eigenfunc- 
tions and the correlation exponent r\ is related to the generalized dimension D 2 of the inverse 
participation ratio rj = d — D 2 . Numerical support for this scaling form has previously 
been obtained from calculations of the two-particle spectral function and the wavepacket 
return probability of real-space and tight-binding models of the quantum Hall effect [|S],|3| as 
well as three-dimensional tight-binding models [|J. 



III. SHAPE OF WAVEPACKETS 

The shape of wavepackets P(r,t) can be obtained by transforming eq. ([!]) back to the 
space and time domain using the scaling form (J|). The long time short distance behavior 
of P(r, t) is thus dominated by multifractal eigenfunction fluctuations characterized by the 
exponent rj. The inverse Laplace transform with respect to u in the limit of q d t/hu ^> 1 is 

1 

P(q,t) ~ (q d t)- D ^ d . (6) 

The shape of the wavepackets is given by the Fourier transform of this expression. The short 
distance behavior, r d <C t/Tip, is thus a power law || 

P(r,t) ~ t - D * /d r D2 ~ d . (7) 

It should be noted that while this is the short distance behavior it can cover many decades 
for long times and hence describes completely the bulk of the wavepacket for long times. 

It also follows from eq. (||) that the return probability P(r = 0,t) decreases as a power 
of time § 

P(r = 0,t) ~t- D2 ' d . (8) 

We notice that the shape of the wavepacket and hence the return probability depends of 
the generalized dimension D 2 reflecting the influence of the multifractal eigenfunction fluc- 
tuations. In contrast to that, the moments of the probability distribution of a wavepacket 



m k{t) = J d d rr k P(r,t) do not depend on D 2 . From the scaling form (||) it follows that 
P(q,t) = p(q d t/hp) and 

m k (t) = J d d rr k P(r,t) = Q d J" drr d ~ 1+k J A^e^P^, t), (9) 

= t k / d Q d fdyy d - l+k J ^e-^p(x), (10) 

where x = q d t/hu and fid is the surface area of the <i-dimensional unit sphere. The moments 
of the wavepacket are integral aspects of the shape of the wavepacket that are not dominated 
by its short distance behavior. Instead due to scale invariance they scale like t k l d . 

The behavior of wavepackets at the mobility edge should be contrasted to the well-known 
metallic behavior of asymptotically Gaussian wavepackets 

-r 2 /(4Dt) 

with constant diffusion coefficient D and moments m^it) ~ t k l 2 . The exponents at the 
critical point differ from these results in that a) the space dimension d is replaced by the 
generalized dimension D 2 due to the multifractal fluctuations and b) the factor 2 is re- 
placed by d reflecting the scale invariance of the conductance that gives rise to factors with 
exponents d — 2 in the critical case. 

It is further instructive to compare our results for the wavepacket dynamics at the mobil- 
ity edge of a disordered system to those obtained for general quantum-mechanical systems. 
Under quite general assumptions it was shown that the return probability decays like t~° 2 , 
where D 2 is the generalized dimension of the spectral measure or the energy dependent local 



density of states |T0] . In the other hand, the exponent of the spatial decay of a wavepacket is 
given by D 2 — d, where D 2 is the generalized dimension of the spatial dependence of the local 
density of states discussed above . Finally, in the absence of multiscaling the moments of 
a wavepacket grow like t k<3 with f3 = D 2 /D 2 0. Our results obtained above are a special 
case of these more general results. They depend only on a single exponent D 2 since at the 
mobility edge the generalized dimensions of the spatial and spectral dependence of the local 
density of states are related to each other, D q = dD q [j3|,|TT| . The origin of this simplification 



is the occurrence of a single relevant length scale L u (eq. (0)), connecting the energy and 
length scales. 

IV. NUMERICAL RESULTS 

In order to numerically study the time evolution of wavepackets for a system defined by 
a Hamiltonian it is necessary to diagonalize at least the part of the spectrum containing 
the support of the wavepacket and to calculate the corresponding eigenstates. Alternatively, 
one can apply the time evolution operator exp(—iHt/h) to a wavepacket. Here the problem 
is that the kinetic and potential energy terms in the Hamiltonian do not commute so that 
a discretization in small timesteps and the use of an exponential decomposition is necessary 



I4IT3II . In any case, the spectrum contains in general extended, critical, and localized states 
at the same time. It is therefore difficult to extract the dynamics at the critical point from 
any superposition of eigenstates of a finite system. Both of these problems can be avoided 



by investigating network models [|i4]-|T6|j. These models are defined by a unitary network 



operator U that can be interpreted as a time evolution operator [|TT|. The disorder as well as 
the energy of the system enters the network operator as parameters. By tuning the energy 
and disorder to their critical values the network operator describes a system at criticality 
with no contributions from localized or extended states. A discrete time evolution in such a 
system is obtained by repeatedly acting with the operator U on a network state ip\ 

m = tfty(o), (12) 

with integer t. Since all eigenstates of U are critical it is possible to start the time evolution 
with a wavepacket maximally peaked at a single site. Again, for a Hamiltonian system this 
is impossible due to the presence of eigenstates with different localization properties in the 
spectrum, limiting the minimum size of a wavepacket. 

In order to investigate the dependence on the dimensionality we study both two- and 
three-dimensional systems. The two-dimensional system is a network model proposed by 
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Chalker and Coddington as a model system exhibiting the integer quantum Hall effect [15 



It has only localized eigenstates except for states at the single energy in the center of a 
Landau level that are critical. The three-dimensional model consist of a stack of coupled 
Chalker-Coddington networks corresponding to a layered system in a strong magnetic field 
[T7J. This system exhibits a true metal- insulator transition with extended and localized 



phases separated by a mobility edge with critical states. 

The behavior of the moments of wavepackets is shown in figs. ([!]) and (|2]). For the 
quantum Hall system in fig. ([!]) the logarithm of mk(t) 2 /k is plotted versus ln£ for k = 
2, 4, 6, 8. Note that in d = 2 the behavior of the moments is the same in the diffusive as well 
as the critical regime. In contrast for the three-dimensional system the exponent is different 
in the critical regime, k/d, and in the diffusive regime, k/2. This is illustrated in fig. ([|), 
where 777.4 {t) is plotted for various energies from the metallic regime (topmost curve) to the 
localized regime (bottom curve) with the critical regime with an exponent of 4/3 in between. 

The decay of the return probability has already been studied in refs. [[|,|J finding agree- 
ment with the expected exponent of D 2 /d. Here we present results for the shape of wave 
packets in a two-dimensional network in fig. (^). The topology of the sample is a torus and 
the time is chosen such that the power law decay is observable up to half of the diameter of 
the system. The observed exponent of 2 — D2 ~ 0.44 agrees well with previous calculations 
of D 2 from other quantities. 



V. CONCLUSIONS 

The time evolution of wave packets at the mobility edge of disordered systems differs 
from the diffusive evolution of wavepackets in good metals. The wavepackets spread for 
all times but their shape never becomes Gaussian. The origin of this behavior is the scale 
invariance of the conductance and the multifractal fluctuations of the local density of states. 
The latter dominate the short distance, long time behavior of wavepackets. The return 
probability decays like t~ D2 / d and the bulk of wavepacket has a power law shape ~ r D 2~ d _ 



7 



The moments of a wavepacket depend only weakly on the multifractal fluctuations in that 
their growth exponent in time, k/d, does not depend on D 2 while the prefactor does. That 
only a single multifractal exponent D 2 occurs and not one for the spatial and one for the 
spectral structure of the local density of states, as is generally the case, can be interpreted 
as consequence of the existence of a single length scale L u connecting energy and length 
scales. 

The numerical study of the time evolution of wavepackets is most economically performed 
for network models compared to models defined by a Hamiltonian. The network models are 
defined by a unitary network operator that provides a natural discrete time evolution. The 
problem then reduces to repeated applications of the network operator. The results of the 
numerical calculations are in agreement with the scaling arguments and with numerical 
calculations based on different methods. 
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FIG. 2. Moments of a wavepacket rrik(t) versus time for a three-dimensional system and differ- 
ent energies showing diffusive (top), critical (center), and localized behavior (bottom). 
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FIG. 3. Probability distribution P(r, t) of wavepackets for a quantum Hall system. The average 
over 5 different realizations of the disorder is shown. The size of the system is 160 x 160 and double 
periodic boundary conditions are imposed. 



11 



